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Abstract. We consider the electron-vibron coupling in suspended nanotube quantum 
dots. Modelling the tube as an elastic medium, we study the possible coupling 
mechanism for exciting the stretching mode in a single-electron-transistor setup. Both 
the forces due to the longitudinal and the transverse fields are included. The effect of 
the longitudinal field is found to be too small to be seen in experiment. In contrast, the 
transverse field which couple to the stretching mode via the bending of the tube can 
in some cases give sizeable Franck-Condon factors. However, the length dependence is 
not compatible with recent experiments [Sapmaz et al. |Cond-mat / 0508270 . 



1. Introduction 

Suspended nanotubes form a interesting and promising system for various nanoelec- 
tromechanical device setups and has been studied both experimentallyjTJ |2l E] and 
theoreticallyjU El El El 13 E]- Because of their large aspect ratio, nanotubes can be 
modeled as simple one- dimensional strings using classical elasticity theory P]. Here we 
study the electromechanical couphng when suspended nanotubes are put in a single- 
electron-transistor (SET) setup. 

When nanotubes are contacted by electrodes they form in most cases contacts with 
a large resistance, which results in Coulomb blockade behavior. This type of single- 
electron-tunnelhng devices have also been fabricated with the nanotubes suspended 
between the two electrodes [HEl Ej- For these devices the interesting possibility occurs 
that a coupling between the electronic degree of freedom and the vibration might show 
up in the current. Such a coupling has indeed been observed in several experiments. 
The first example LeRoy et al.[2j observed phonon sidebands with both absorbtion 
and emission peaks, which were taken as evidence for the radial breathing mode being 
strongly excited and thus behaving essentially as a classical external time-dependent 
potential. 

In the quantum regime the electron-vibron coupling leads to steps in the IV 
characteristic, similar to the well-known Franck-Condon principle. This has been 
seen in a number of single-molecule devices [Tn| [TT| and is well understood in terms 
of rate equations[T21 113 HH HE]- Recently, similar physics was observed in suspended 
nanotubes [3] where the vibrational energy suggested that the sidebands were due to 
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Figure 1. Illustration of suspended nanotube device. 



coupling to a longitudinal stretching mode. However, the coupling mechanism was 
unclear. It was suggested in Ref. [Sj, that the electric field parallel to the tube coupled to 
the nuclear motion. Here we will argue that due to screening in the tube the longitudinal 
electric field is too weak to excite the longitudinal mode, and instead point to the non- 
linear coupling between the traverse and longitudinal mode as the possible coupling 
mechanism. 

The paper is organized as follows. First, we discuss in section |21 the electrostatics 
of the charged suspended nanotube, followed by an account in section El of the elastic 
properties of the hanging tube. In section IHl the modes of the tube are quantized 
and the Franck-Condon overlap factors are calculated. Finally, conclusions are given in 
section [3 

2. Electrostatics of charge nanotube 

We now discuss the electrostatics of the charge suspended nanotube. First we consider 
the longitudinal electric field, and then the radial field. 

2.1. Electric field parallel to the nanotube 

For a metallic tube there would, of course, be no electric field in the longitudinal 
direction. However, the tube has a finite screening length due to the kinetic energy. 
We have analyzed this situation using the following density functional for the total 
energy of a nanotube with linear dispersion 

Ti 1 
F[p] = — / dxvp{x)[p{x)]'^ + - / dx dx' p{x)V{x,x')p{x'), (1) 
2 Jo 2 Jo Jq 

where vp is the Fermi velocity and V{x,x') is the effective ID potential for a cylindric 
conductor. Details about the interaction and the solution are given in Appendix A. One 
should include screening due to both gate and source-drain electrodes. The gate can be 
included as a metallic plane at some distance h, and the source-drain electrodes as a 
section of the wire with = 0. See figure 2 for an illustration of this. 

To minimize the energy functional in equation ([Q) under the constraint that a charge 
of one electron is added to the tube, we add a Lagrange multiplier 

Fi[p, \] = E[p] + \ (^j^ dxp{x)-e)j. (2) 
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Figure 2. The model used to find tlie electric. The electrodes are represented by a 
ID lead with vp — 0. 



First, we minimize with respect to p and find 
5F 

hvpp{x) + / dx'V{x,x')p{x')+X = 0, (3) 



6p[x] 

and then with respect to A: 

-Q^ = J dxp{x)-e = 0. (4) 

These two hnear equations are readily solved numerically. Once the solution is found, 
the electric field is given by 



d 



L 



eEJx) = —— / dx' V(x, x')p(x'). (5) 
ox Jo 

The important parameters in this solution are the aspect ratio of the tube, i.e., ^, and 
the strength of the interaction 

^ = T-^- (6) 

Aireohv-p 

For typical parameters, one has r = 10 — 20, while the aspect ratio is 200 — 2000. The 
distance to the gate is not important, as long as it longer than the length over which 
the electric decays, which is typically the case. Our numerical solution gives an electric 
field comparable to the results of Guinea ^H] (see also reference llT,), which results in 
an electric force of order eE^^ ~ 10~^ N, for typical nanotube devices. However, this 
field is limited to a small region near the contacts, and therefore the total effect of it is 
small. 

For r ^ 1, we can in fact make a simple approximation which quit accurately 
describes the full numerical solution. The electric field can related to the charge density 
by differentiating the condition Q with respect to x: 

d 

eE^{x) = h— [p{x)vf{x)] . (7) 

Because p{x) changes little along the wire, we may set p ^ 1/L and we thus obtain 

eE^{x)^^[6{x)-6{x-L)]. (8) 

The width of the delta function will be given by microscopic details of the interface 
between the tube and the contact, i.e. a length scale of the order of the radius of the 
tube itself. This length scale we denote by xq. 
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2.1.1. The electrostatic force in the longitudinal direction The term in the Hamiltonian 
describing the interaction between the electron charge density and the nuclear system 
is 

-f^el-vib.x = - y" dxdx p{x)V{x,x)pn{x'), (9) 

where Pn{x) is the density of the positive ions in the tube. The force per length acting 
on the mechanical degrees of freedom is therefore given by eE^pn. In terms of the 
displacement field defined below in equation (fTHjl . Hc\-vih,x becomes 

-f^ei-vib,x = -po j dx p{x)V{x,x')[dx'u{x')\ = -epQ j dxE^{x)u{x). (10) 
2.2. Electric field perpendicular to the nanotuhe 

To find the electric field in the radial direction we model the nanotube as a distributed 
capacitor similarly to Sapmaz et al. [IHj. We include capacitances to the electrodes Ci 
and to the gate 

C = Cl+Cg, (11) 

where the capacitance to the gate is 

with c being the distributed capacitance of a tube over a plane. 

To find the total charge on the tube, we write the total electrostatic energy as 

W=^- gA$/e, (13) 

where C is total capacitance with charge q and A$ is the difference between the 
nanotube and the electrode work functions. (Here we neglect the effect of the source, 
drain and gate voltages, because they are considerably smaller than A0.) The optimum 
charge is thus 

go = rioe = A<l> C/e. (14) 

For single-walled carbon nanotubes the workfunction is about 4.7 eV |19t 1^ I21j. while 
for gold it is 5.1 eV. For typical devices C ~ 10"^'' F and hence uq ~ 30. The electrostatic 
energy is used in the following section to calculate force acting on the tube. 

2.2.1. The electrostatic force in the transverse direction Below we solve for the 
distortions of the wire due to the electrostatic forces. The force in the direction 
perpendicular to the charged wire (denoted the z-direction) is given by 

dW 5c„ 



^9 



where 



dCg 5z{x 

5c, 



z=0 



Cg \ e'^nl 



-a 



8z{x) 



C ) 47reo/iL2' 
2'nen Cl 



_ 

T r 1 -1 / 

z=0 



h[cosh-\h/R)]' 2neohL^' 



(15) 



withCg = Lc(0). (16) 
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Figure 3. The coordinate system used to describe the hanging tube. A point on 
the tube, which before the distortion was at (a;',0) has after the deformation the 
coordinates [£,[x'), z{x')). 



3. Elastic properties 

In this section, we discuss in detail the elastic properties of a suspended nanotube. Most 
devices are made by under-etching after depositum of the nanotube, which is done at 
room temperatures. Therefore, since the thermal expansion coefficient of nanotubes 
is negative (100 nm wire expands a few nm from room temperature to 4 K) it seems 
unlikely that there is a large tension in tube unless it becomes heavily charged [THj. 
When the tube is charged it is attracted towards the metallic gates and leads. 

The radial force, off course, couples to the breathing-type modes, which, however, 
have too large energies (~ 30 meV) to be relevant for the low-voltage experiment in 
reference jS]- Here we are interested in the lower part of the excitation spectrum and 
disregard this mode. We are left with the bending and stretching modes. The energy 
of the bending is typically much lower than those of stretching modes[7], and therefore 
we treat these as purely classical, which means that we will solve for the bending mode 
separately and then via an anharmonic term this solution acts as a force term on the 
longitudinal mode. 

We thus consider two possible mechanism for coupling to the stretching mode: 
either directly via the longitudinal electric field discussed in section 12.11 or through the 
perpendicular field, section which bends the tubes and hence also stretches it. 

3.1. Elasticity theory of a hanging string 

Assuming the tube to lie in a single plane and thus having no torsion in the longitudinal 
direction, we can describe the distortion of the bend tube by (^(x), z{x)), where x G [0, L] 
runs along the unbend tube in absence of external forces, see figure |H1 (If the tube has 
some slack, x is a curve linear coordinate along the equilibrium position the hanging 
tube.) This means that a point along the tube, which before was at (x, 0) after the 
deformation has the coordinates {S,{x) z{x)) . The total elastic energy of the tube is 
then follows from standard elasticity theory of an elastic string [22] 




(17) 
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where C, is the hnear strain of extension, R is the radius of curvature, A = irR"^ the area 
and / = 7tR^/4 the moment of inertia, E is Young's modulus, and k\\, k±_ the external 
forces, and where we have defined the longitudinal displacement field as 

u{x) = C,{x) — X. (18) 

The linear extension of an infinitesimal element between x and x + dx is 



Cix)dx = ^/{^{x + dx) - ^{x)y + {z{x + dx) - z{x)y - dx, (19) 



or 



C(x) = (^y/[l+u'{x)]' + [z'{x)]^ - l) . (20) 
The linear extension elastic energy is thus 

W^lin=— / dx[^[l+u'{x)f + [z'{x)f - l) . (21) 



The curvature contribution is determined in a similar way. First, the unit tangential 
vector is 

{i\x),z\x)) 



(22) 



(23) 



which then gives the square of the radius of curvature as 

and then the curvature contribution to the elastic energy finally becomes 

_El f\ iz'{x)u"ix)-il + u'{x))z"ix))' 
2 io {[l + u'{x)]^+[z'{x)Y)^ ■ 

3.2. Weak distortions 

Since we are interested in small deflections, we expand the two elastic energy expressions 
for small z and u. For VTun, we obtain to third order in u and z 

dx([u'{x)f +'^^^^^+u'{x)[z'{x)f]. (25) 



lin o , , L- v"/j ■ 4 

Here we note that the last term couples the bending and stretching modes. For the 
curvature contribution, we find to the same order 

ET ET 

P^cur = — / dxil- Au'ix)) [z"ix)f ^— dx [z"ix)f . (26) 
Jo Jo 

Again, there is a term which couples the two modes. However, for nanotubes this term is 

much smaller than the last term in equation (j25j) . because it smaller by a factor (R/L)^, 

and hence we have neglected the coupling term in PFcur- 
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4. Solution for the bending of the tube 

As mentioned above, the bending mode itself has a resonance frequency too low to 
be seen tunneUing spectroscopy (~ 100 MHz [7]) (even when under tension due to the 
charging of the wire ) , which means that it can be treated as a classical degree of freedom 
and Franck-Condon type physics is not involved. In contrast, the longitudinal stretching 
mode has been seen in the single-electron-transistor (SET) device [3.. and here we wish 
to calculate the Franck-Condon coupling constants for this mode. Therefore we take 
the following approach: first we solve for the bending mode classically and then insert 
this as an external force acting on the longitudinal mode. The differential equation for 
z{x) is 

lEz"" - ^{z'fz" = k. (27) 

This equation cannot be solved analytically. One approach is to approximate {z'Y by 
the average value, which is corresponds to assuming constant tension in the wire |23j. 

Below we solve for the bending function z{x) in two regimes: the linear and the 
non- linear regime. Once we know z(x), we will be interested in the change of z{x) due 
to tunnelling of a single electron. For large uq, the relevant change is thus 

ze{x) = (28) 
auo 

This change will then couple to the longitudinal mode via the coupling term in 
equation (j^ . 

4.I. Linear regime 

For weak forces we can simply neglect the non-linear term in equation ()27|1 . and with 
boundary conditions z{0) = z{L) = z'{0) = z'{L) = the solution is 



24E/ V LJ \L 

The shift in zq{x) due to the charge of a single electron is according to equation (j28|) 
given by 

^o.w^^^(4r(^r^l-l)^m^ (30) 



12neohEA\Rj \C J \ LJ \L. 

For a tube with R = 0.6 nm, L = 1 fim, and E ^ 10^^ Pa and a device with h = 200 
nm, no = 50, and Cg/C = 0.5, the maximum distortion is of order a few nm. 

The linear approximation is valid when the second term in equation ()27j) is much 
smaller than the first. Using the solution in equation (I^Uj). this translates to the condition 

10-«m7^)%<l. (31) 



^RJ \EI ^ 

For the typical parameters used below, the number on the left hand side of (|HT|) is < 1 
and therefore the linear approximation is only marginally valid. 



Electron-vibron coupling in suspended nanotuhes 



8 



4-2. Non-linear regime. 

For larger distortions the non-linear term in equation p7jl becomes important. In the 
strongly non-linear regime, we can neglect the first term and we have 



AE 



z'fz" 



AE d 



'\3 



2 ' " ' ~ 6 dx 
with boundary condition z{0) = z{L) = 0. The solution of this equation is 



z[{x) 





r 


X 1 


[ea, 




L~2 



1/3 



Sign ( ^ 



X 



(32) 



(33) 



In this non-linear regime, the change in the slope the bending function zi{x) due to a 
single electron charge, is 

/n \ 2/3 / .2 



C ) \2T:tQhAnQ 
5. Distortion of the longitudinal mode 



1/3 



X 

L 



1/3 



Sign 



X 



(34) 



Due to the electrostatic forces the equilibrium position of longitudinal displacement field 
u{x) shifts. Since the forces are small the displacement is going to be small. For the 
tunnelling overlap factors, the important point is, however, whether these displacements 
are large compared to the quantum uncertainty length, which later is seen to be of order 
of pm. In this section, we calculate the classical displacements of u{x). One example is 
shown in figure IH 



5.1. Distortion due to the longitudinal electrostatic force 

The displacement of the longitudinal mode is readily found from its equation of motion. 
The displacement due to the longitudinal electric field follows from 

EAu'^ix) = k\\{x) = -eE^Po, (35) 

with boundary conditions u{0) = u{L) = 0. With forces concentrated near the contacts 
as in equation (jH)), there is an abrupt change of u{x) at x = and x = L. See figure |21 
red dashed curve. 



5.2. Distortion due to the transverse electrostatic force 

Once we have solved for z{x) in equation ()28|). we can find the force that acts on 
the longitudinal diplacement field u by inserting the solution into the last term of 
equation (j2SI), and then identify the force k± in equation (fTTj) . This gives 

The size of the displacement follows from the balance between this force and the strain: 
EAu'^ix) = k^{x) ^ u'^{x) = z'{x)z"{x), (37) 
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Figure 4. The solutions for the shifted longitudinal mode u{x) due to the parallel 
(red dashed curve) and perpendicular (blue curve) electric fields. We have used typical 
parameters as in section [6.31 



which together with the boundary condition, Mo(0) = uq{L) = 0, gives the solution 



«o(x) = -- dy [z\y)r + — dy [z'iy)]'. (38) 

One example is shown in figure |21 (blue curve) . 

In the next section, we analyze the Franck-Condon overlap factor due to this 
electrostatic distortion of the stretching mode. 

6. Quantum mechanics and Franck-Condon overlap factors 

The longitudinal eigenmodes of a nanotube modeled as a ID elastic medium follows 
from the Hamiltonian 

H= r dx( ^ + ^{d,,u{x)A , (39) 



JO \ '^Prn 2 / 

where pm is mass density per unit length, and u and p are conjugated variables, 
i.e., [u{x),p{x'y\ = ih6{x — x'). In order to diagonalize H, we introduce the Fourier 
transforms: 

\ /Trnx\ ^ ^ \/2 . /TTnx\ , 

n=l 

/Tcnx\ ^ ^ r 1 /Tmx\ ^ 
p(x) = — 2^ sm (^-— j p„, pn = V2 dzsmi^-—jp{x), (41) 

n=l "^O 

where and Pn obey [-UrajPn] = i^- Now H transforms to 

Pi , AE 



^ = EU^ + Tr(™)'^n ' M = p^L, (42) 

n=l ^ ^ 

The Hamiltonian ()42p is easily diagonalized by 

K = {an + ai) , Pn = j^\J\{0'i-d^) . (43) 
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where and are usual annihilation and creation operators, and 



, AE -n AE -KVs „ h , , ,^ 



With these operators, the Hamiltonian ()42|) becomes 

(45) 



n=l ^ ^ 



As we saw in the previous sections, additional terms in the Hamiltonian appears 
due to the force generated by the tunnelling electron. These are included next. 

6.1. Coupling due to longitudinal electric field 

The longitudinal electric field E^ gives rise to a coupling Hamiltonian (see equation (fTUj)): 

(46) 



Hei~vih,\\ = -epo / dx u{x)E^{x) = V'un/n 

■^0 n=l 



where 

-L 



/„,, = -ep„V-2l dx sm {^) E^.) « -(2.v/2)!^5^iMS, (47) 

for n even and zero for n odd. 

6.2. Coupling due to the capacitative force 

The transverse force leads to the following term in the Hamiltonian (see equation ()25|1 ) : 

EA °° 

^cl-vib,± = ^ / dx U'{x)[z'^{x)]^ = y2'^nfn,±- (48) 
^ n=l 

where 

fn,l. = dx COS (^-^j [<(x)]l (49) 

6.3. Franck- Condon overlap factors 

The tunnelling of an electron leads to a displacement of the equilibrium displacement 
field according to equations pUjl . fPHj) . Each mode represented by Un is shifted by the 
amounts 

C,a = ^|f^, a = {l±). (50) 

This allows us to calculate the Franck-Condon parameters, which for each mode n 
expresses the overlaps between the new eigenstates around the new equilibrium positions 
and the old one. This parameter is defined asflH] 

9n,a = \i^-^^' ■ (51) 



2 



'0,n 
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The size of g determines the character of the vibron sidebands in the /^^-characteristics, 
such that for g there are no sidebands, for g of order one clear sidebands are seen, 
while for ^ 1 a gap appears in IV characteristic [T^ ITH] . 

For the Franck-Condon parameter due to the parallel electric field, we thus have 

'-'^^ = ^^[-AEr) ■ (''^ 

Using a;o ~ 1 nm, L = 500 nm, fp = 10® ms~^, R = 0.6 nm, and E = 10^^ Pa, we find 
92,]] ~ 10~^. This is clearly too small a coupling constant to explain the experimental 
findings in reference [O]. 

The coupling constant due to the perpendicular electric can be expressed explicitly, 
for the case of small z{x), using equation ()29|1 . The integral in equation ()49|1 can then 
be performed and we obtain 

MQL^ (kjr\ ^ (nV-^40)' 

and zero for n odd. Using as defined in equation (fT5|l . and typical parameter for 
single wall carbon nanotube devices: R = 0.6 nm, L = 1 /im, E = 10^^ Pa, ~ 10 — 200 
nm, uq = Lc(O)A0 ~ 30, Cg/C = 0.1 — 0.75, we find the maximum gn factor to occur for 
n = 4. However, for this range of parameters we also find (?4 -C 1, unless the geometry 
is such that 

^°"'>0.1. (54) 



9n,± = — z — ( ) — o 13 14 — ' ^ e've^' (53) 



Even though we can get a significant coupling, the condition ()54j) does not seem to be 
compatible with the experimental realizations in reference Even more so because 
the coupling strongly strength, equation (j^HI), depends strongly on the length of wire, 
which is not seen experimentally. 



7. Conclusion and discussion 



We have considered the electromechanics of suspended nanotube single-electron- 
transistor devices. When the charge on the tube is changed by one electron the 
resulting electric field will distort the tube in both the longitudinal and transverse 
directions, and both these distortions couple to the stretching mode. We have calculated 
they consequences for the coupling constant for vibron-assisted tunnelling expressed 
as the Franck-Condon factor. This is expressed in terms of the ratio of the classical 
displacement to the quantum uncertainty length. Even though both are in the range of 
picometers, the effective coupling parameters, g, turn out to be small for most devices. 

Because the screening of the longitudinal electric field is very effective, the dominant 
interaction seems to be the coupling via the bending mode. However, only if the tube 
is very close to the gate do we get a sizeable ^f-parameter. This could indicate that in 
the experiment of Sapmaz et al.0| the suspended nanotube has a considerable slack, 
which would diminish the distance to the gate. Further experiments and more precise 
modelling of actual geometries should be able to resolve these issues. 
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Appendix A. Electric field in a charged nanowire 

Appendix A.l. The effective ID interaction 

Here we derive in more detail the result show in section El The interaction between two 
charges on points x and x' is 

^^|0(rx)n0(rl)p 



Vix-x ) = dr^ dr\ , (A.l) 

where 0(rx) is the wavefunction is the perpendicular direction. By modelling the tube 
as a cylinder, 0(r_L) = 5{r — R)/2nR, we get 

f^^ de 

v/(x-x')2 + (2i?sin(^/2)2 

k(^-^^) , (A.2, 



V{x - x') = — - / 
27reo Jo 



27r2eo^(x-x)2 + 4/?2 ~ x'^ + AR\ 

where K is the complete elliptic integral of the first kind. With screening due to a gate 
at distance h ^ R, the interaction is 

Vscrix - x') = V(x - x') . (A.3) 
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